Spin- Velocity Correlations of Optically Pumped Atoms 
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We present efficient theoretical tools for describing the optical pumping of atoms by light prop- 
agating at arbitrary directions with respect to an external magnetic field, at buffer-gas pressures 
that are small enough for velocity-selective optical pumping (VSOP) but large enough to cause sub- 
stantial coUisional relaxation of the velocities and the spin. These are the conditions for the sodium 
atoms at an ahitude of about 100 km that are used as guidestars for adaptive optics in modern 
ground-based telescopes to correct for aberrations due to atmospheric turbulence. We use spin and 
velocity relaxation modes to describe the distribution of atoms in spin space (including both pop- 
ulations and coherences) and velocity space. Cusp kernels are used to describe velocity-changing 
collisions. Optical pumping operators are represented as a sum of poles in the complex velocity 
plane. Signals simulated with these methods are in excellent agreement with previous experiments 
and with preliminary experiments in our laboratory. 



I. INTRODUCTION 

In this paper we discuss efficient ways to model opti- 
cally pumped atoms in a regime where velocity-selective 
optical pumping (VSOP) is possible, but where coUi- 
sional rates with buffer gases are too high to permit the 
use of models for cooling and trapping in the near ab- 
sence of collisions. This is the regime of sodium guidestar 
atoms. These naturally-occurring layers of sodium atoms 
at altitudes of 90-100 km above the Earth's surface [ij 
are illuminated by ground-based lasers, and the return- 
ing photons are used to measure the relative retardation 
of wave fronts across an optical aperture. This retarda- 
tion information can be used with a deformable mirror 
to correct for the aberrations from atmospheric turbu- 
lence and to allow the receiving optics to produce a more 
nearly diffraction-limited image of astronomical objects. 

The performance of guidestar systems is limited by 
the loss of atoms from the most strongly backscatter- 
ing spin sublevels and velocity groups. The most im- 
portant reasons for these losses are: collisions with the 
residual atmospheric gases, which transfer atoms from 
strongly absorbing to weakly absorbing spin sublevels, 
or which shift the atoms into velocity groups that are 
not in resonance with the pumping light; Larmor preces- 
sion of the spins away from strongly absorbing orienta- 
tions if the geomagnetic field is not parallel to the direc- 
tion of the laser beam; and unwanted optical pumping 
into weakly absorbing sublevels. The powerful model- 
ing methods discussed here make it easier to explore the 
parameter space of these processes and to optimize the 
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performance of guidestar systems. These methods also 
provide a more realistic and numerically convenient way 
to model laboratory experiments with VSOP of atoms in 
low-pressure buffer gases. In contrast to previous work 
on this topic, for example references , we account for 
the full hyperfine structure of real alkali-metal atoms, we 
show how to use spin-relaxation modes [1, Q to incorpo- 
rate into the model the complicated spin relaxation of Na 
guidestar atoms due to collisions with paramagnetic oxy- 
gen atoms, and we use recently developed cusp kernels 
yf\ to realistically and efficiently model velocity-changing 
collisions. 

This article is organized as follows. In Section [III we 
introduce the Liouville space of the coupled spin and ve- 
locity distributions of the atoms. In Section Hill we intro- 
duce spin-relaxation modes Sj d, Q to describe the spin 
distributions, and we introduce the concept of conjugate 
spin-mode indices, j and j*. The amplitudes \xj) of the 
spin modes in velocity space provide a complete descrip- 
tion of the spin and velocity polarization of the atoms. 
The little-known Liouville conjugate operation, denoted 
by the superscript and the transposition operator T for 
Liouville-space operators are discussed in Section IIIIBI 
Using Liouville conjugates reduces the numerical com- 
puting requirements by nearly a factor of 2. 

In Section IIVI we show how to describe velocity dis- 
tributions with velocity-relaxation modes |w„) 8, 9]. We 
use the velocity modes to show that a simple transforma- 
tion of the widely used Keilson-Storer kernels [13] leads 
to much more realistic and useful cusp kernels 0] for 
describing velocity-changing collisions. In Section |V] we 
present a simple model for the transition from collision- 
free, ballistic flight to any container walls at very low 
buffer-gas pressure to diffusional wall losses at higher 
pressure. We sketch how the relative sizes of the laser 
beam and the cell affect these processes. In an Appendix 
we show how to deduce the rate of velocity-changing col- 
lisions, 7vd from the spatial diffusion coefficient D and 
the smallest-nonzero eigenvalue ai of the collision oper- 
ator Avd with the little-known formula (|A.1[) . In ((5^ 
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of Section I VII we show that spin-changing and velocity- 
changing cohisions cause the spin- mode amphtudes \xj) 
to relax exponentially in time at the rate Kj, where Kj 
is a kernel in velocity space. 

In Section [VIIl we introduce a velocity-dependent op- 
tical pumping operator ^op, which we write as the sum 
of poles in the complex-velocity plane at locations de- 
termined by the laser frequency and the optical Bohr 
frequencies. The pole expansion facilitates velocity aver- 
ages in terms of Faddeeva functions (Voigt profiles) [lli • 
The poles have "residue matrices" that are independent 
of the laser frequency and atomic velocity. In Section 
VIII, we show that the steady-state mode amplitudes \xj) 
generated by the combined effects of optical pumping, 
spin relaxation, and velocity relaxation can be written in 
terms of Green's functions, G,- oc K~^. We show that 
if the kernel Kj is a cusp kernel or linear combination 
thereof, Ky^ is also cusp kernel or linear combination 
thereof. Being able to invert cusp kernels in closed form 
greatly simplifies the numerical evaluation of the mode 
amplitudes \xj)- This simplification is not possible with 
Keilson-Storer kernels or any other collision kernel that 
we know of. We present an explicit, series solution (|134p 
for the mode amplitudes \xj) in powers of the pumping 
light intensity, with particular emphasis on the first-order 
solution (|TT7l) . 

Finally, in Section IX, we use (|147p to demonstrate 
how the methods we present are in excellent agreement 
with existing VSOP experiments. We also describe a new 
type of magnetic-depolarization experiment that can be 
carried out under laboratory conditions and readily inter- 
preted with the powerful modeling methods described in 
this paper. Such experiments would provide much more 
detailed experimental information about the nature of 
velocity-changing collisions. 



Here x{^) denotes a square g^^^ x g^^^ matrix m 
Schrodinger spin space for ground-state atoms, where the 
dimension of the spin space is g^^^ — {2S + 1)(2/ -I- 1). 
The electronic-spin quantum number of the ground-state 
atom is S, and the nuclear-spin quantum number is /. 
The total probability for the atom to have some velocity 
and be in some spin state must be unity, so we must have 

Tr / xix)dx = 1. (4) 



.{9} 



A. Energies and energy basis states 

It will be convenient to describe the atoms in terms 
of the energy eigenstates |^). These are defined by the 
time-independent Schrodinger equation 



i7{f>|/.)=i?^|M), 



(5) 



which determines the energy shifts i?^ of the basis states 
l/i) from their center of gravity due to hyperfine in- 
teractions and externally applied magnetic fields. The 
spin Hamiltonian H^^^ for the ^Si/2 ground state of an 
alkali-metal atom is traceless and includes hyperfine cou- 
plings of the nuclear and electronic spins to each other as 
well as their couplings to an externally applied magnetic 
field. The energy sublevels |/2) and energy shifts Efi of 
optically-excited atoms are given in like manner by 



(6) 



The Bohr frequency for an optical transition from the 
sublevel to the sublevel |/2) is 



Efi — 



(7) 



II. THE DENSITY MATRIX 

For optical pumping at low buffer-gas pressure we need 
to account for both the spin-polarization and the veloc- 
ity V of atoms along the light beam. We introduce a 
dimensionless velocity. 



V 
VD ' 



(1) 



where the most probable speed vd along the laser beam 
is given by 



Here lo^^^^ = ck^'^^^ — 2nc/X^'^^^ is the mean value of the 
frequencies averaged over all possible combinations 
of the sublevel labels, /2 and /i. The corresponding spatial 
frequency is k^*^^^ and wavelength is A''-'^^'^ . 

For the low geomagnetic fields of interest to us, it will 
sometimes be convenient to use low-field labels of the en- 
ergy sublevels /i ^ f^ra^. Here, / denotes the approx- 
imate total spin angular momentum quantum number 
of the sublevel and m is the exact azimuthal quantum 
number along a quantization axis defined by the external 
magnetic field. 



Vd 



2kBT 
M ■ 



(2) 



Here T is the absolute temperature, fc^ is Boltzmann's 
constant, and M is the mass of the atom. 

We write the incremental density matrix dp for the 
spin-polarized, ground-state atoms with velocities be- 
tween X and X + dx as 



dp = x{x) dx. 



(3) 



B. Liouville space 

We use a generalization of the Liouville-space for- 
malism of the recent book Optically Pumped Atoms by 
Happer, Jau and Walker [I4I (which we will refer to 
as OPA) for handling the large amount of information 
needed to describe the spin- velocity correlations of opti- 
cally pumped atoms. One of the best early descriptions 
of Liouville space is given in the book by Ernst et al. [l^ . 
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To describe the (dimensionless) velocity x, we can use 
Tlx evenly spaced sample velocities, 



c„^ , with Xk+i - Xk = Sx. (8) 



Then the velocity-dependent spin polarization can be de- 
fined by the elements of the spin density matrix x of ([21) i 
which we will call "spin-velocity correlations," 



(9) 



We will think of {xk\Xi-iiy) ^-s the projection onto the 
velocity-space basis vector \xk) of the abstract, velocity- 
space column vector |xaii^)- It will be convenient to rep- 
resent the total density matrix for spin-velocity space as 
the abstract, Kronecker-product column vector 



(10) 



where l^iy) is the Liouville-space representation of the 
spin density matrix basis element We turn now 

to the time evolution of |$). 



where the Bohr frequencies of the ground-state atoms are 

(14) 



In (|12p spin-changing collisions occur at a rate 7sd and 
the details of the spin relaxation are described by the 
dimensionless matrix operator 



\^lv){^lv\Asd\^J.'l''){^^'v'^ 



(15) 



The specific form of the damping operator Agd for var- 
ious collisional interactions is discussed in Chapter 10 
of OPA [l2|. Regardless of the particular details, for 
the evolution described by ([TT|) to conserve the number 
of atoms, the spin-evolution operators must satisfy the 
constraints 



((so|r = ((so|i/© = ((so|Asd-0. 



(16) 



The equilibrium left eigenvector, about which we will 
have more to say below, is 



(17) 



III. SPIN-DAMPING 

There is negligible correlation between spin relaxation 
and velocity relaxation for laboratory VSOP experi- 
ments, where all the spin relaxation is due to collisions 
of polarized atoms with the cell walls, or for sodium 
guidestar experiments where the spin relaxation is is 
almost all due to binary spin-exchange collisions with 
oxygen molecules, and where the electron spin, but not 
the nuclear spin of the Na atom may flip. We there- 
fore take the spin-relaxation processes to be independent 
of velocity-relaxation processes and we write the rate of 
change of ([T0| due to the hyperfine interaction, the exter- 
nally applied magnetic held, and spin-changing collisions 
as 



d_ 
dt 



<i>)=.-^r|H®lxM.) 



(11) 



Here the spin-damping operator is independent of veloc- 
ity and is given by 



-H 



© 



- 7sd^ 



sd- 



(12) 



The evolution due to internal hyperflne couplings of the 
electron and nuclear spins as well as their interactions 
with externally applied magnetic fields is given by the 
Liouville-space Hamiltonian, a "commutator superoper- 
ator" as described by (4.85) of OPA 



(13) 



A. Spin Modes 

In this section we discuss how to handle the compli- 
cated spin relaxation of sodium guidestar atoms due to 
gas-phase collisions with O2 molecules and O atoms. For 
the weakly-relaxing buffer gases used in laboratory ex- 
periments, spin- relaxation collisions in the gas are slow 
enough to be neglected, and the spin relaxation is almost 
entirely due to collisions with the walls. An especially 
convenient basis for the spin polarization under condi- 
tions of strong collisional relaxation in the gas is provided 
by the spin-relaxation modes. To our knowledge, spin 
modes were first introduced by M. A. Bouchiat [H, |B| to 
discuss the curious multi-exponential decays observed in 
velocity-independent optical pumping. The spin modes 
are the right eigenvectors of the spin-damping operator 

(da, 



r|«j) = 7j|si). 



(18) 



Here the symbol \sj) denotes a g^^^^ x 1 column vector 
in Liouville space, formed from a g^^^ x g^^^ matrix Sj 
of Schrodinger space, as described by (1.2) and (1.3) of 
OPA [l2| , by placing each column of sj below the one on 
its left. We will enumerate the modes such that 



l7o|-0<|7i|<|72|<...<|7«.-i|. 



(19) 



^1/ 



For degenerate spin modes \sj) ^ \sj') ^ \sj")... with 
7j = jji = "fj"... one is free to take linear combinations 
of spin modes with other distinguishing properties, often 
the angular momentum. A simple example of longitudi- 
nal spin modes is shown in Fig. 10.5 of OPA [12||. 
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In the absence of any pumping mechanisms, the spin 
density matrix of the atoms will relax to the thermal 
equilibrium state 



.0 = \e 



The partition function Z is 



Z = Tr 



(21) 



The evolution operator F must cause no changes to the 
steady state, or 



r|so) = o. 



(22) 



The left eigenvectors of F satisfy an eigenvalue equa- 
tion analogous to (|18l) with the same set of eigenvalues, 



((sjl7j- 



(23) 



We will use the symbol ((sj | to distinguish a left eigenvec- 
tor, the solution of (|23l) . from the Hermitian conjugate 
(sjl — |sj)^ of the right eigenvector |sj), the solution 
of (|18p with the same eigenvalue 7j. It is necessary to 
make this distinction because there may be right and left 
eigenvectors for which ((sj | is not a simple multiple of 
(sjl — \sj)'^ . The right eigenvectors |sj) are analogous 
to the three primitive vectors of a crystal. For crys- 
tals of low symmetry, which correspond to non-normal 
operators, the primitive vectors can be non-orthogonal. 
The left eigenvectors ((sj| are analogous to the recipro- 
cal primitive vectors. Except for rare, singular combi- 
nations of the parameters of (|12l) , the right eigenvectors 
\sj) span the spin space of the atom. As long as the right 
eigenvectors are linearly independent we can choose left 
eigenvectors such that 



((sj|sfc) 



and ^ I 



I {4 



(24) 



Using the left and right eigenvectors we can write the net 
spin-damping operator as 



i 



(25) 



High-temperature limit. For our applications, the 
thermal energy ksT will normally be so large compared 
to the energy differences of ground-state sublevels, that 
it will be an excellent approximation to simplify the equi- 
librium spin mode (PU)) to 



iM |iM) 
So = -j^ or \so) - 
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l{9} 



(26) 



where the unit operator for the Schrodinger ground state 

is i^«> = EIm)(m|. 



B. Hermitian conjugates of spin modes 

According to (4.4.2) of OPA [13], F must be Liouvil- 
lian, 



F^ = F, 



(27) 



to keep the density matrix Hermitian as it evolves in 
time. From this requirement, we can derive some identi- 
ties that dramatically reduce the amount of computation 
necessary to simulate a guidestar or VSOP scenario, and 
that simplify the form of the evolution equations. The 
Liouville conjugate F^ of F is defined by 

pt ^ Tr*T. (28) 
The transposition operator of (|28l) can be written as 

T = J2Wfi){fi,.\. (29) 

Substituting ^ into (JHl) and using ^ we find 

TT*T\sj)=jj\s,). (30) 

Multiplying ^ on the left by T, noting that = 1, and 
taking the complex conjugate of the resulting equation we 
find 



TTls* 



7*T|4). 



(31) 



Here we have defined the complex-conjugate column vec- 
tor by 

i^„.\s*) = {^,\s*W) = ms,W)]* = iiHsjW- (32) 

From pT|) we see that 



(33) 



is an eigenvector of F with eigenvalue 7* . With some care 
in the case of degeneracies, where several spin modes have 
the same eigenvalue 7^ , we can define the modes such that 



r|s,)-l4) or l4)-k,.) 
with 7* = 7j*. 



(34) 



Since taking two Hermitian conjugates restores the orig- 
inal matrix, we must have 
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J- 



(35) 



Using the expressions above we find that an alternate 
expression for the transposition operator of (|29l) is 



(36) 



Hopefully, the use of the superscript * to denote the index 
j* of the mode sj* = s], will not be confused with the 
symbol for complex conjugation of a number. The indices 
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for the Us = {21 + 1)(2S'+ 1) modes are the real integers, 
j = 0, 1, 2, . . . , ris — 1 and the mode indices j* are the 
same real numbers in some permuted order. Both Sj 
and its Hermitian conjugate s| are square matrices in 
the ground-state subspace of Schroedinger space; \sj) and 
|sj) are column vectors in Liouville space, but |sj)^ = {sj\ 
is a row vector in Liouville space with elements {sjlfiiy) — 
{fiiy\sj)*. 

The matrix elements of the Liouville conjugate 
TM*T of a spin-space matrix M are 

(^i/|M*|kA) = {i^^\M\\k)*; 
is,\Mi\sk)^(s,.\M\sk.r. (37) 

C. Spin-mode expansions 

For numerical work, it will be convenient to replace the 
spin- velocity basis vectors |/LtJ/) (^\xk) and ($j {xk \ by 
basis vectors with a spin-mode part and a velocity part, 

\sj)<S)\xk) and (sj\<g){xk\- (38) 

We can expand the density matrix |<i>) of (fTO|) . on the 
basis vectors ([38]) to find, 

m^Y.\'^j)^\xj)- (39) 

The coupled and uncoupled velocity amplitudes of ((39)) 
and (jlOp are related by 

For Hermitian density matrices ^''^ = $ we must have 
\x*) = \Xj'), (41) 

or 

{xk\x*)^{xk\Xjr ^{xklXj')- (42) 

The elements {xk\xj) of population modes with j* = j 
must be real. 



IV. VELOCITY-DAMPING 

We assume that gas-phase collisions cause the density 
matrix PH)) to change in velocity space at the rate 

d 

^|*) = -7vd5Ikj)®^vd|Xj)- (43) 

j 

Here 7vd is a characteristic velocity-damping rate, which 
can be deduced with the aid of (jA.l[) from the spatial dif- 
fusion coefficient D and the smallest, non-zero eigenvalue 
ai of the collision operator Avd- The collision operator 



is often written in terms of a velocity-changing collision 
kernel Wvd as 

Ad = - W^vd. (44) 

The velocity-space unit operator is 1^"^^ = ^ \xk){xk\- 

Atom conservation implies that that the velocity- 
evolution operators of (j44p must satisfy a constraint anal- 
ogous to (|16p for the spin-evolution operators, which we 
write as 

(volA^d = 0, (45) 
ivolW^d = (vol (46) 

The equilibrium left eigenvector ((uqI is analogous to ((so| 
of P?]) . and is given by 

ivo\=J2('^k\. (47) 

k 

The equilibrium right eigenvector \va), conjugate to the 
left eigenvector {vo\ of (|47)) . is analogous to |so) of (|20)l . 
and is given by the Maxwellian distribution 

M-^5]K)e-^^. (48) 

The analogs of the equilibrium mode constraint ()22|) are 

A.d\vo) = 0, (49) 
Wvd bo) = bo). (50) 

A. Velocity modes 

In analogy to ([T8|) and ([23| . we assume that the 
velocity-damping operator Avd has a spectrum of right 
and left eigenvectors \vn) and ((t;„| corresponding to the 
eigenvalue a„ such that 

Wvdbn) = n7„|i;„) and ((w„|Wvd = ((^^n|^7„, (51) 

Adbri) = Q;„|t)„) and ((t)„|Avd = (I'nlan, (52) 

with 

n7„ = 1 - «„. (53) 

One of the first clear examples of velocity modes was 
given by Snider Q. There will be the same number 
Ux of independent eigenvectors \vn) as there are ve- 
locity sample points Xk- The eigenvalues are real and 
non-negative. They can be numbered by the integers 
n = 0, 1, 2, . . . , rij; — 1 such that 

zuq = 1 > TUi > W2 > • . • , and 

ao = < ai < a2 < (54) 

In analogy to ([24]) . we assume that the left and right 
eigenvectors can be chosen to be orthonormal and com- 
plete, so that 

WW) = (Jnfe, and E \Vn)(Vn\ = l^"'^- (55) 
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The values of the summation index are n = 0,l,2,...In 
analogy to we write 



Wvd = ^ tJ7„ 1 1;„ ) {vn I , and A^d ^^an\Vn) (Vn 



(56) 



Keilson-Storer kernels 



A convenient model kernel Wa = Wvd for the velocity- 
changing collision kernels Wvd of ([44|l was introduced by 
Kcilson-Storer (KS) [10]. In the KS model, it is assumed 
that a group of atoms, all having the the same sample 
velocity Xk along the laser beam, are transformed by an 
ensemble of single collisions, with various impact param- 
eters and orbital planes, into the distribution of final ve- 
locities x.. 



Wa{Xj,Xk) = -^{Xj\Wa\Xk) 



(57) 



This newly formed Gaussian distribution is centered at 
the dimensionless velocity axk- The "memory" parame- 
ter a and width b are 



< a < 1 and b = \/ \ - a} . 



(58) 



Snider [8| has shown that the KS kernel can be written 
as the eigenvalue expansion [§| 



W^a = ^a'>„)((i;„|. 



(59) 



The amplitudes of the left and right eigenvectors can be 
chosen to be 

„ I X Hn{xk) J / I N SxHnixk)e~'''i 

iVn\Xk) = —==- and [Xk\Vn) ^ 



V2"n!7r 



(60) 

Here Hn{x) denotes the nth Hermite polynomial. The 
KS eigenvectors are independent of the memory param- 
eter a. Then the KS eigenvalues of ([5^ are 



and 



1 



(61) 



With their simple analytic form KS kernels ([57)1 are 
convenient for numerical work, and they clearly satisfy 
the normalization constraint (j46p and the Maxwellian 
constraint (|50p . However, KS kernels with a single mem- 
ory parameter a do not give very good approximations to 
kernels inferred from experimental observations, for ex- 
ample, those of Gibble and Gallagher [TJ] or to kernels 
modeled from realistic interatomic potentials, like those 
of Ho and Chu [H. 

Real collisions occur for a large range of impact pa- 
rameters, or as large numbers of partial waves in a quan- 
tum treatment of the scattering. "Head-on" collisions 
with small impact parameters will produce large changes 
in velocity and can be approximately modeled by KS 



kernels that are close to the strong-collision limit, with 
a = 0. "Grazing- incidence" collisions with large impact 
parameters will be much more frequent but will produce 
small changes in velocity. They are better modeled by 
KS kernels with awl. 

With these facts in mind it has often been proposed 
that a superposition of KS kernels would be a better 
model, but suggested superpositions have been somewhat 
inconvenient for numerical modeling. McGuyer et al. Q 
have introduced the "cusp kernel," a special superposi- 
tion of KS kernels that is even more convenient for mod- 
eling than KS kernels, since cusp kernels and their su- 
perpositions can be readily inverted to find steady-state 
velocity distributions. It is more difficult to invert KS 
kernels or other collision kernels that have been used in 
the past. Cusp kernels and their superpositions are also 
more similar to kernels inferred from experiment (l^ or 
from realistic interatomic potentials [l5l |. 



C. Cusp kernels 

Cusp kernels describe an ensemble of collisions, with 
each sample collision described with a KS kernel of mem- 
ory parameter a. The probability to find the the memory 
parameter between a and a + da is Ps{a)da, where the 
probability density is 



Ps{a) — sa 



s-l 



(62) 



We will call the parameter s the "sharpness," and the 
corresponding velocity-changing collision kernel of (|44|) 
will be denoted by W^d = Cs- The sharpness s will 
be a real, positive number for velocity-changing collision 
kernels Wvd, but for "resolvent kernels" W , which we will 
discuss below, the sharpness can have an imaginary part. 
An explicit expression for the cusp kernel is 



Cs 



Jo 



s\Vn)iVn\ 



The eigenvalues of ([53|) are 



Wn = and 

s + n 

In the limit n — > oo we have 



n 



s + n 



Woo = and aoo = 1, 



(63) 



(64) 



(65) 



for both cusp kernels and KS kernels (|5T|) . 

Morgan and Happer [9| have shown that the matrix 
elements of (1631) can be summed to give 

Cs{xj,xu) = ^^^e"'='i?,(-a:<)i?,(a;>). (66) 

Here T{s) is the Euler gamma function, Xy is the greater 
of the two variables Xj and Xfc, and a;< is the lesser. The 
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"right function" Rs{z) can be represented with the power 
series convergent for all finite z, 



^{-zT 



^^n!2-«r(i + f-§)- 



(67) 



For sufficiently large sharpness, \s\ ^ f , one can evaluate 
the cusp kernel (j66p with the asymptotic expression 



2\nCs{xj,Xk) = xl-x]-2\xj-xk\^/2's + \n(^^ . (68) 

A comparison of KS kernels and cusp kernels is shown in 
Fig.ffl 
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FIG. 1. (Color online) Keilson-Storer kernels Wa of (f57|l with 
representative memory parameters a, compared to cusp ker- 
nels Cs of (|66p with sharpness parameters s that have the 
same expectation value, (a) = J^^ Ps{a)ada = s/(s + 1), of the 
memory parameter. As a — >■ and s — !■ 0, the KS and cusp 
kernels generate Maxwellian distributions. 

As shown by McGuyer et al. [7[, one can fit collision 
kernels inferred from experimental measurements very 
well with the superposition of a few (m « 3) cusp kernels: 



vd 



fc+1 



fkCs^, with 



fc=l 



fk 



I. 



(69) 



Each cusp of sharpness Sfc makes a fractional contribution 
fk > to the overall damping operator Wvd • An example 
of a multicusp collision kernel from McGuyer et al. [7j 
is shown as the curve labeled p = in Fig. ^Bji.. This is 
the kernel that was used for the representative models of 
experimental data in Figs. [5l |6] and |8l The other parts 
of Fig. [5] are discussed in Section IVIII Al where we show 
how to find multicusp resolvent kernels from multicusp 
collision kernels. 
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FIG. 2. (Color online) Resolvent kernels W versus buffer-gas 
pressure p for the initial velocity y — 1.37. Here, the collision 
kernel Wvd of (|69p is a three-cusp fit from Ref. 0] to a mea- 
sured kernel for Rb in He from Ref. [l^]. The superposition 
parameters are fj = [0.13,0.37,0.50] and Sj = [ 7 8, 27.2 500]. 
(a) As described by the expressions of Section IVIII Al with 
increasing buffer-gas pressure the resolvent kernel W evolves 
from a low-pressure, single-collision limit, W = Wvd, to 
a high-pressure limit, which is a Maxwellian distribution 
W — \vo){vo\. (b) With increasing pressure the weights gk 
of the resolvent kernel tend to the limit gk = [1, 0, 0]. (c) The 
sharpnesses of the resolvent kernel diminish with pressure, 
and in the high pressure limit, ri 0, give a Maxwellian dis- 
tribution for the only resolvent cusp with appreciable weight, 
gi ~ 1. (d) The pressure-dependent rates 7vd of (|A.1|I and 
7w of (|79p are for a — 1 mm, 6 = 5 mm, T = 60 °C, and 
Do = 0.15 cm^s~^ at 760 Torr. The transition from low- 
pressure to high- pressure conditions occurs when jw ~ 7vd. 
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V. WALL-DAMPING 



B. Cylindrical cells 



Most existing data on velocity-selective optical pump- 
ing has come from laboratory experiments, where the 
atoms either experience no buffer gas collisions at all, or 
they collide with buffer gases like He, Ne, Ar, Kr, Xe, 
N2 , or H2 , where the rate 7sd of spin-changing collisions 
in ()12|) is orders of magnitude smaller than the rate 7vd 
of velocity changing collisions in Under these con- 

ditions, the damping rates 7^ of the spin modes can be 
well approximated by 7^ — for population modes with 
j* — j, and by jj = iujj for coherence modes with j* j 



and Bohr frequency u>j 



-UJj 



Under these labora- 



tory conditions, the effects of the buffer gas on velocity- 
selective optical pumping are characterized by the rate 
7vd of velocity-changing collisions and by the effective 
rate 7w of collisions with the wall. Walls with no spe- 
cial coatings destroy the spin polarization of impinging 
atoms and release unpolarized atoms with a Mawellian 
distribution of velocities. 

For analyzing the physics of Na guidestar atoms 
there is no need to consider wall collisions, but the O2 
molecules and O atoms of the buffer gas at an altitude 
of 90-100 km have rates 7sd for spin-changing collisions 
that are comparable to or larger than the rates 7vd for 
velocity-changing collisions. Under these conditions, the 
damping rates 7-,- of the spin modes have relatively large 
real parts due to spin-changing collisions. 



A. Walls 

In analogy to ([TT]) and (|43|) we write wall relaxation as 
d 



dt 



(70) 



The effective collision rate of atoms with the walls, 7w, is 
a real, positive number. We assume that every atom that 
hits a wall sticks, and is replaced by a different atom that 
evaporates with no spin polarization and with a Maxwell 
distribution of velocities. The wall depolarization opera- 
tor is simply 



— 1 



(71) 



Here we use (|17p with (|47p to write the equilibrium spin- 
velocity row vector for the full spin-velocity space as 

i0\^iso\®iv^\. (72) 
The corresponding equilibrium column vector. 



I0) = |so) ® \vq), 



(73) 



is the density matrix for atoms with no spin polarization 
and a Maxwell distribution of velocities. For the modi- 
fications needed to account for walls that only partially 
depolarize the spins and lead to shifts of the coherence 
frequencies, see the work of Wu et al. [l6| . 



As a semiquantitative example of how to estimate wall 
relaxation rates 7w, let us suppose that a pump laser fills 
a cylindrical volume of radius o on the axis of a cylindrical 
cell of radius b > a. Let x be some spin-mode amplitude 
of the density matrix that is pumped at a rate S inside 
the pump laser beam. We assume axial symmetry so 
that X — x('") depends only on the distance r from the 
cylinder axis. Then the diffusion equation is 



d_ 

dt 



X-DV\ 



for 
for 



r < a 

a < r < b 



(74) 



The spatial diffusion coefficient can be measured experi- 
mentally and the results are often expressed as 



(75) 



where p is the gas pressure and Dq is the diffusion coeffi- 
cient at the reference pressure poj which is usually 1 atm. 

One can integrate the steady-state version of (TMl) with 
the boundary condition x(&) = and with x and dx/dr 
continuous at r = a to find the solution 



X 

s 



(a^ - r^) 
AD 



iii„(^ 

2D \ a 



if r < a, 



X 

s 



2D \ r 



if a < r < b. 



(76) 



(77) 



For pressures high enough for the diffusion equation to 
be valid, we define the mean spin polarization x sampled 
by the probe beam, and the mean time for an atom 
to be lost from the probe beam by 



1 



X 



x{r)2'Krdr = St^ 



(78) 



Substituting (|76l) into ([78| and adding a representative 
free flight time, a/vo, for the atom to escape the pump 
beam in the limit of very low pressures, we find 



1 



a 

VD 



a 

8D 



1 



41n( - 

a 



(79) 



VI. EVOLUTION IN THE DARK 



Summing the rates of change of the density matrix 
(pg)) . from wall collisions (f70|) . from hyperfine interactions 
and spin-relaxing collisions (jlip . with F given by p5|) . 
and from velocity relaxing collisions (j43p , with A^^ given 
by ((56|) . we find that the evolution rate from all sources 
except optical pumping is 



d 



^l*) = -E(^w + 7,)|sj)®IXj) 



+ 7vd I Sj ) «) ^vd I ) + 7w 1 0) (f 1 Xo ) 



(80) 
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For laboratory and guidestar experiments with weak 
pumping light, the unpolarized part of the density matrix 
will have a Maxwellian distribution 



-1 



+1 



+2 



Ixo) = \vo). 



(81) 



For spin-polarized parts of the density matrix with j ^ 0, 
we can muhiply on the left by ((sj| ® to find 



d_ 
di 



\Xj) = -(7w + Ij + 7vdAvd)|x,) = -K,\xj). (82) 



Here we have introduced a damping kernel for the jth 
spin mode, 



Kj = ^ljn\Vn)iVn\- 



(83) 



The kernel Kj includes the evolution of the spin due to 
gas-phase and wall collisions, hyperfine interactions, and 
precession in an external magnetic field. Also included 
in Kj are velocity-changing collisions. The characteristic 
relaxation rate for the nth velocity mode of the spin mode 
j is 



Ijr. 



7w +7j + a«7vd- 



(84) 



Here 7w is the effective collision rate of polarized atoms 
with the wall, 7j is the sum of relaxation due to gas- 
phase collisions (= for most laboratory experiments) 
plus a factor iujj for the Bohr frequency ujj of the spin 
mode |sj), and Q;„7vd is the relaxation rate of the velocity 
mode \vn) due to gas-phase collisions at the rate 7vd- 
According to (|34|) and (|84|) the coherence damping rates 
can be complex, but they must satisfy the identity 

l*n = l3'n. (85) 

A special case of (|84|) is 700 = 0. 



VII. OPTICAL PUMPING 

The basic processes involved in optical pumping are 
sketched in Fig. [3l We consider monochromatic pump- 
ing light of temporal frequency lo = 2t:u and spatial fre- 
quency k — Lu/c, propagating along the unit vector ^. 
We assume that the light intensity is large enough to 
cause substantial spin polarization of the ground-state 
atoms, but that it is not so intense that it produces 
substantial population of the excited state. Most lab- 
oratory experiments on velocity-selective optical pump- 
ing, and most guidestar systems are in this regime. It 
is a straightforward extension to account for the satu- 
ration of the optical pumping of the ground state from 
intense, repetitively-pulsed lasers. For laboratory experi- 
ments, the atoms will also be subject to a weak, counter- 
propagating probe beam, usually a small fraction of light 
from the source of the pump beam. As the laser fre- 
quency changes, resonant changes in the attenuation of 
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FIG. 3. Energy sublevels of an alkali-metal atom, resonant 
frequencies vjj, and representative transitions for absorption 
of light (depopulation pumping), or spontaneous emission of 
light (repopulation pumping). In this example the pumping 
light is linearly polarized parallel to a small, externally applied 
magnetic field. Because of the Doppler shift, the resonance 
frequencies depend on the atomic velocity and on the direction 
of propagation of the light. 



the probe beam occur for laser frequencies where the ve- 
locities of atoms spin-polarized by the pump beam also 
have the right Doppler shift to resonantly absorb light 
from the retro-reflected probe beam. The classical elec- 
tric field of a monochromatic light beam for an atom 
located at the position r at time t is 



E 



^ —i{<7kC,-r—ujt) 



(86) 



The direction index a — ±1 will be taken to be ct = 1 for 
the pumping beam. For a probe beam that propagates in 
the same direction as the pumping beam we will also have 
(7 = 1, and for for a counter-propagating probe beam we 
will have a — ~1. The amplitude E of the probe beam 
may have a different polarization that than the pump 
beam. The optical intensity / (in units of ergcm~^s~^) 



is given in terms of the field amplitude E of 



by 



27r' 



lEl' 



(87) 



As discussed in connection with (5.52) of OPA ^T^], the 
interaction of the light with the atoms can be well ap- 
proximated by the sum of conjugate interaction matrices 
V + V''', where 



V 



-Dt • E and V'' 



D E*. 



(88) 



The component of the electric dipole moment opera- 
tor D, when operating on a ground state basis j/u), pro- 
duces an excited-state basis \p,), and vice versa for D. As 
shown in connection with (6.7) pf OPA [T^j, the unsatu- 
rated optical coherence generated between excited-state 
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and ground-state sublevels of the atom is proportional to 
the matrix, 



W = V.jE 



{eg} 



(89) 



In ([89]) we use the "dot-slash" symbol (./) to denote 
element-by-element division of the matrices with the 
same dimensions. 

Resonant velocities. From (5.84), (5.88), and (5.92) 
of OPA 12] we see that the energy-denominator matrix 
that occurs in is 

^p.'lf^ = fi(CTfc^^^^'yfc-w+Wp^)-ifi-7c = fi-7D (o-a^fc - Zfifj.) ■ 

(90) 

The kth discretized component of the atomic velocity 
along the pump beam is Vk ^ v ■ — XkVn- The photon 
momentum is very nearly ahk^'^^^^, where the spatial 
frequency is k^'^3\ The complex resonant velocities are 
determined by the optical frequency uj and by the optical 
Bohr frequencies of the atom, and given by 



(91) 



The real and imaginary parts of the resonant velocities 
are 



Id 



and Uc 



— , (92) 

7-D 



with assume that the damping rates 

associated with the optical coherences /I/i are the same 
and given by 



1 



7c 



2r{«} 



+ 7o 



(93) 



where r^*^^ is the natural radiative lifetime of the ex- 
cited atom. The collisional damping rate 70c of the op- 
tical coherence is nearly negligible for most VSOP sit- 
uations. In modeling, the parameter 70c will be used 
to approximately account for the frequency linewidth of 
the laser and for slight misalignment of the pump and 
retro-reflected probe beam. 

The pumping will be resonantly enhanced for veloci- 
ties Xk that are close to the resonant velocities Xnn , that 



is, for Xk 



Similarly, the probe absorption will be 



resonantly enhanced for velocities as close as possible to 
the resonant velocities aXjy^, that is, for x^ ~ oxp^. If 
the probe beam propagates parallel to the pump beam 
(cr = 1), every resonant velocity for the pump beam will 
also be resonant for the probe beam. However, if the 
probe beam is counter-propagating with a = —1, the 
resonant conditions can be satisfied for the same velocity 
group only if the laser detuning is such that two different 
optical transitions, /2/i and vv, have equal and opposite 
resonant velocities. 



■-MM 



(94) 



This is the condition for saturated-absorption resonances 
with counter-propagating laser beams. In Fig. |3]we show 



how the resonant velocities and their negatives —Xuv 
depend on laser detuning for a '^^K atom in a magnetic 
field of 1 G. 

Residue matrices. Using (|90p we can write (|89p as 



aXk - Zn 



(95) 



For each optical transition between a ground-state sub- 
level fi and and excited-state sublevel /I, we have defined 



a gi"^ X g 
ment 



{g} "residue matrix" with a single nonzero ele- 



S--S ^'^'^ 



D 



(96) 



The residue matrices are independent of the laser fre- 
quency w and the atomic velocity x^. They give the 
relative contributions to the pumping rate of the optical 
transitions from sublevel \^) to 



A. Depopulation pumping. 

A natural rate to characterize velocity-dependent op- 
tical pumping is the maximum possible pumping rate. 



7mx 



(97) 



for a hypothetical isotope with no externally applied 
magnetic field, with vanishing hyperfine coupling coeffi- 
cients, with no spin polarization, and with zero velocity. 
Such hypothetical atoms would have only a single optical 
Bohr frequency = to^^^^ , and their absorption cross 
section for light of frequency uj would be independent of 
optical polarization and given by 



7c 



where the maximum cross section is 



(98) 



(99) 



Here f^^e} jg i^j^g oscillator strength, is the classical 
electron radius, c is the speed of light, and the optical 
coherence damping rate was given by (l93l) . 

We will write the evolution of the mode amplitude \xj) 
due to depopulation pumping as 



dt 



(a^felXj) = -lm^^{S]\AAp{xk)\si){xk\xi)- (100) 



We have neglected the small changes in velocity Xk 
caused by absorption and emission of light. The depop- 
ulation pumping operator Adp can be written as 



^dp(a;fe) = A^pixk) + ^dp(^fc) 



(101) 
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Density matrix amplitude (arb. units) 




Laser detuning v — i/^^f} (GHz) 

FIG. 4. (Color online) Resonant velocities ±a;p^ of (|92p for atoms moving along a pump beam (the solid lines sloping up to 
the right) and for a retroreflected and attenuated probe beam (the dotted lines sloping downward to the right) as a function 
of the detuning of the laser frequency from the center of the 770 nm Dl hue of ^^K atoms in a magnetic field of B = 1 G. 
Both the pump and probe light are linearly polarized parallel to the magnetic field, as indicated in Fig. [S] The absorption 
frequencies from Fig. [3] to which each group of resonant velocities correspond are denoted by vj^. For a group of unpolarized 
atoms with dimensionless velocity x = 0.75, we display the absorption cross sections a^(y,x) of a pump beam (the solid line), 
and probe beam a-{v,x) (the dashed line). Plotted as soUd lines along the vertical axis at the center figure are the first-order 
population shifts {x\x\}^^ of H148P as a function of velocity x for the ground-state sublevel with low-field quantum numbers 
/ = 2 and m = 2, produced by a pump beam with a frequency detuning of GHz. Plotted as a dashed line on the same scale 
is the velocity dependence of the factor 5R((so|^dp(3;)|/U.M) of (|107[) for the probe beam. Collisions in a 10 mTorr buffer gas at 
T = 50 °C have transferred atoms in the sublevel from their four resonant excitation velocities to other velocities, some of 
which are in resonance with the probe beam. For these curves, a = 0.1 cm, 6 = 0.5 cm, and Do = 0.1 cm'^s"^. The dot-dash 
line is a Maxwellian distribution of velocities, (2;|uo) of (|48() . Strong VSOP resonances may be observed at laser frequencies 
that make pump and probe resonant velocities equal, within the envelope of the Maxwellian velocity distribution. 



The Liouvillc-conjugate matrix is defined in terms 



of Adp by (Hg). We can use (6.12) of OPA [ij 
to find 



R. 



and dMl) 



(102) 



B. Absorption rates and optical cross sections 

The rate of depletion of ground-state atoms by depop- 
ulation pumping is 



(105) 



where the residue matrix is given by 



-ptyylw} 



(103) 



X 9 



{g} 



As discussed in (4.4.1) of OPA [H], a g^^} 
Schrodinger-space matrix M is transformed into into a 
corresponding g^^^^ x g^^^^ "fiat matrix" of Liouville 
space by the Kronecker product 



where (ST) is the rate at which spin-polarized atoms ab- 
sorb or scatter light for a light beam of intensity /, and 
the corresponding absorption cross section of the spin 
polarized atoms is (a). Combining (llOSp with (jlOOp we 
find 



(ST) = {5f) + {Sty. (106) 



where 



(^f) = 7mx^((so|idp(Xfc)|s,)(Xfc|Xj). (107) 
jk 



= l^ff} (g) M. 



(104) 



Mean pumping rate. We define the mean optical pump- 
ing rate, 7op of the atoms as the value of the pumping 
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rate (ST) for unpolarized atoms with a Maxwellian dis- 
tribution of velocities, and therefore with the spin mode 
amphtudes 



(a^felXj) = {xk\vo)Sjo, 



(108) 



Then we can use (fT08)) with (ITOTl) and ((T02]) to write the 
amphtude of the mean pumping rate as 



"Top — 7mx ^ ^ 



= *V^7mx^((so|i?l7^|so)u;(zpp). (109) 



The integral over the velocity distribution (that is, the 
sum on k in (|109l) ) can be written in terms of the Fad- 
deeva function [17], which for Q{z) > is given by 



w{z) 



dt 



t- z 



(110) 



The Faddeeva function is a superposition of Lorentzians 
with a Gaussian distribution of resonance frequencies. It 
is often called a Voigt profile, and it can be evaluated 
with a very efficient computer algorithm due to Weide- 
man 11]. 

The mean pumping rate 7op of (|109p depends on the 
laser frequency through the parameter z^^ and is pro- 
portional to the laser intensity because of the factor 7nix- 
For large magnetic fields, 7op also depends on the laser 
polarization. The equilibrium absorption cross section 
for unpolarized atoms with a Maxwellian distribution of 
velocities to absorb light is 



7op 



(111) 



where the light intensity / was given by (I87p . 

Repopulation pumping. In analogy to (|100[) , the rate 
of change of the spin-mode amplitudes due to repopula- 
tion pumping is 



^(a^felXi) = 7i 



nix 



(4.91) of OPA [T^l, only the diagonal elements of the o- 
dot transform X® of a Schrodinger-space matrix X are 
nonzero, 

(Aii^]X®]/ii/) = (mIXIj/), and 
i^^l,\X'^~'\^,ly) = {fi\X\iy)-\ (115) 

The spontaneous emission matrix is given by (5.50) of 
OPA nl as 



(116) 



The sum extends over the projections Dk = • D of 
the dipole moment operator, where x/j is a unit vector 
along the fcth Cartesian axis of a spatial coordinate sys- 
tem. From (5.45) of OPA [12] we find that for pumping 
through the first excited ^Pj state, the squared ampli- 
tude of the dipole operator that occurs in the denomina- 
tor of ()116|) is given in terms of the spontaneous, radiative 
decay lifetime t^'^^ of the excited atom by 



w 



2 _ hc^i2J+l) 



(117) 



Net optical pumping. The sum of the depopulation 
pumping (llOOp and repopulation pumping (jll2l) is the 
net optical pumping, 

d 

Ol^^klXj) = -lmx'^isj\Aop{xk)\si){xk\xi)- (118) 



The amplitude of the net optical pumping operator Aop 
is 



A-op{xk) = ^ ' ■ 



axk - zn 



with R^i:^R^;-R^^^. 

_ (119) 
As shown in (6.72) of OPA [IJ] the net optical pumping 
operator satisfies the constraint that no atoms are created 
or destroyed by optical pumping, 

((soliop(a;.) = or (sol^^ ^ 0. (120) 



In analogy to (|102p the amplitude of the repopulation 
pumping operator is 



A-pixk) = 



^rp 

axk - Zn 



(113) 



Using (6.71) and (6.35) of OPA [H, we find 



(114) 

Here H^''^ excited-state Hamiltonian, and t^''^ is the 
radiative lifetime of the excited atoms. As defined by 



VIII. STEADY-STATE SOLUTION 



Adding the contribution ([82)) for relaxation in the dark 
to the contribution (jll8l) from net optical pumping we 
find 

d 

OllXj) = -KjlXj) - lmK'Y\Xr)isj\Aop{Xr)\si){Xr\xl)- 

Ir 

(121) 

For j 7^ 0, the steady-state solution of (|12ip is 

\Xj) = ~PjYGj\xk)isj\Aop(xk)\si)(xk\xi)- (122) 

Ik 
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Here we have introduced the dimensionless Green's func- 
tion 



Gj = ^jooK, \ 



(123) 



and a dimensionless optical-pumping parameter for the 
jth spin mode, 



Pj 



7m 
7jc 



(124) 



proportional to the light intensity. The characteristic 
damping rate 7^00 = 7w + 7j + 7vd was given by dUl) . 
The parameter pj and the amplitude \xj) decrease as 
|7joo| increases. For hyperfine coherences, the imaginary 
parts of 7jooj the hyperfine Bohr frequencies, are so large 
compared to |7jc3o| for non- hyperfine coherences, that it 
is a good approximation to neglect the hyperfine coher- 
ences entirely and retain only Zeeman coherences and 
population modes. 



(I127p . From a partial-fraction expansion of D{n)/Nj{n) 
we find that the weights of the resolvent cusps are 



9jk 



lim 



l]ool]o{n + rjk)Din) 
IvdTjkNjin) 



(130) 



One can show that the fractional weights gjk sum to 
unity. 



E 

k=l 



gjk 



(131) 



An example of the pressure dependence of multicusp re- 
solvent kernels is shown in Fig. [21 In the low-pressure 
limit, with p — ^ 0, the resolvent and collision kernels 
coincide, that is, gj^ and rj^ — )■ s^,. In the high- 

pressure limit with p — > 00 all of the weight is transferred 
to the least sharp cusp, and the sharpness of this cusp 
approaches zero with increasing pressure. 



A. Green's functions for multi-cusp collision 
kernels 

For the multicusp collision kernel of (|69p the damping 
rates (l84l) become 



where jjo = 7w + 7j- The denominator polynomial is 
simply 

D{n) ^ {n + si){n + S2) ■ ■ ■ {n + s,„). (126) 



7jO 



1 

oc -. 



l]ooY.kUk/sk) P 



(132) 



A cusp kernel with sharpness r^i <^ 1 produces a 
Maxwellian distribution from any initial distribution of 
velocities. 

Complex Green's functions. Since 7* = 7^* , we see 
the complex conjugate of the Green's function is 



(133) 



where j* is the conjugate index to j, defined by (l34|) . 



B. Series solution 



Different spin modes may have different numerator poly- 
nomials 



^ji'f^) = ljaD{n) + 7vd ^ 



fknPjn) 
n + Sk 



= 7ioo(« + rji){n + rj2) ■ ■ ■ {n + rjm), (127) 

where the damping rate 7^00 = 7w + 7j + 7vd is the limit 
of dHll) for n ^ 00. Using with §^ and (fT^ we 

find that the multi-cusp Green's function is 



7jO 

The multi-cusp resolvent kernel is 

^J=J293kCr, 



(128) 



(129) 



Here C^^^ denotes a cusp kernel ([55)1 of sharpness rjk, 
given as one of the roots of the numerator polynomial 



For non-equilibrium modes with j ^ we can write the 
mode amplitudes as a power series in the light-intensity 
parameter of (|124p . 



\xj) = Ep; 



(134) 



The optical pumping parameter pj for the jth spin mode 
was given by (|124p . Special cases of ()134p that follow 
from ([5T|) are 



\xf) = S,o\vo), and 1x^0 = ^«o|t'o). (135) 
Substituting ([TM]) into (IT^ we find for j > and n > 0, 



J2Gj\xk)isj\Aopixk)\si) 



(136) 



First- order spin polarization. For the first-order spin 
mode, we can use ()136|) with ()135p to find 



(137) 
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where 



1^1^^) = - ^Gj\xk)isj\Aop{xk)\sQ){xk\vQ) 



(138) 



We have used (I128P to write the amphtude as the sum 
of a part coming from atoms that have had only wah 
cohisions, 



i^k\xrJ) = -j2 



is,\Ki;\so){xk\vo) 



Xk - Zn 



(139) 



and a background or "pedestal" from atoms that have 
had gas-phase collisions, 



ix,\xy/) = i^5](x.|W^,>.)(a:.|x]:0- (140) 
7jO V 



Physical insight can be gained by considering the first- 
order population shifts 



IX 



{1}^ - 



(141) 



Noting from (gB]) that ivo\Wj = ((uo|, we multiply 
on the left by X^iC^oki) to find 

Hlxjg^) = — ((«o|x£^). (142) 

According to (I142p . the collisional background area 

(wolxjg^) and coUision-free area ((wolxjw^) the ratio 

of the velocity-damping rate due to gas-phase collisions, 
7vd, to the gas-free spin damping rate due to wall colli- 
sions, 7jo = 7w + 7j- 

Absorption rate. Using (jl34p we can write the absorp- 
tion rate (|107p as a power series in the light intensity. 



(<5r)=^(5rW), 



(143) 



where 



((5r^">) = ((5f^">) + ((5f^">)*, (144) 
and the amplitude of the nth-order absorption rate is 
(^f{">) =7,„,^p7((so|idp(xfc)|s,)(,-- 



XklXj (145) 



jk 



The zeroth-order rate is 



^ 7op, 



(146) 



with the amplitude 7op of the mean pumping rate given 
by p09)) . The first-order rate is 



(5f{l>)=^2nHi((so|idp(Xfe)|s,)C 



Xk\X 



(147) 



For experiments where only population imbalances 
and no coherences are created by optical pumping and 
7sd ^ 7vd , we can write (|14ip as a sum of contributions 
from each ground-state energy sublevel — |/to), 

(jf = —iso\A^Axk)W){xk\x]!^). (148) 

For computing the absorption rate of a probe beam in 
a pump-probe experiment, /probe should be used in ((97)) 
for computing 7mx, and the resulting (ijr'f^^) must be 
multiplied by the ratio /pump/Vobo- 

Representative examples of the factors 

iso\Adpixk)\^J-^J.) and {xk\xi]}) (|148p are shown 
in Fig. Ill 



IX. COMPARISON WITH EXPERIMENT 

The main purpose of this paper is to describe efficient 
ways to model the optical pumping of Na atoms under 
conditions similar to those of Na guidestar atoms. It has 
proven difficult to simulate conditions of Na guidestar 
atoms in laboratory experiments. It is relatively easy to 
carry out experiments at the same buffer-gas pressures as 
those experienced by guidestar atoms with non-reactive 
buffer gases like He, Ne, Ar, Kr, Xe and N2. Unfortu- 
nately, the spin-damping rates of these gases are many 
orders of magnitude smaller than the velocity- damping 
rates. Because residual air at the ^ 100 km altitude of 
the Na layer is still approximately 20% O2 molecules by 
volume and can even contain a few percent of dissociated 
O atoms, these paramagnetic species will cause the spin- 
damping and velocity-damping rates to be of comparable 
magnitude. Although Na atoms can have many binary 
collisions with O2 molecules in the upper atmosphere, 
with negligible probability for a chemical reaction, the 
walls of a laboratory container quickly catalyze oxida- 
tion of Na atoms by O2 gas. But the following examples 
show that the modeling methods work very well with 
laboratory VSOP experiments. 



A. Modulated Circular Dichroism of Na 

A particularly convenient way to investigate collisional 
effects on velocity-selective optical pumping was intro- 
duced by Aminoff et al. [l^, who measured the effects 
of Ne buffer gas on the saturated absorption resonances 
of Na atoms. The Na in their experiments was held in 
a cylindrical glass cell and spin polarized by resonant, 
589.8 nm Dl laser light, pumping along a small, 1 G 
magnetic field with alternating circular polarization. In 
the limit of weak pumping light, this pumping scheme 
produces an alternating orientation (polarization of mul- 
tipole index L — 1). The polarization is detected as 
the modulated attenuation of a weak, counterpropagat- 
ing probe beam of fixed circular polarization. We will call 
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this a modulated circular dichroism (MCD) experiment. 
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FIG. 5. (Color online) A MCD signal modeled with (IT48t fo r 
measurements like those shown in Fig. 4 of Aminoff et al. 
The modeling parameters were temperature T = 150 °C, 
optical-coherence linewidth 7oc/27r = 10 MHz, longitudinal 
field i? = 1 G, bufi'er-gas pressure p = 30 mTorr, spatial 
diff'usion coeflicient Do ~ 0.6 cm^s~^, laser beam diameter 
a — 3.5 mm, and a cell diameter fe = 10 mm. The cusp-kernel 
parameters of Fig. |2]were used with no adjustment. 
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FIG. 6. (Color online) Magnetic depolarization signal pre- 
dicted with p47p for a transverse field of B = 0.1 G, but with 
all other modeling parameters the same as those of Fig. El The 
"collisional background" of the saturated absorption signals 
is greatly modified and even reversed in sign. Such experi- 
ments could be used to determine the optimal cups kernels 
for describing velocity-changing collisions. 

In analyzing their MCD signals, Aminoff et al. [l8| 
found that KS kernels gave unsatisfactory fits to ob- 
servations, and that much better fits could be obtained 
by adding a phenomenological, narrow Gaussian kernel, 
given by their Eq. (32), to a broad KS kernel. The re- 
sulting two-term kernel does not produce a Maxwellian 



distribution of velocities in thermal equilibrium, that is, 
it does not satisfy the fundamental constraint ([50)1 . The 
two-term kernel also cannot be inverted conveniently to 
obtain a Green's function, and in general it required the 
evaluation of a "cumbersome double sum" [l^ . 

In Fig. [S] we show the difference in absorption rates 
p47p of the probe beam for right- and left-circularly- 
polarized pumping light. Of the 64 possible spin modes 
|sj), only two have amplitudes \xj) that differ apprecia- 
bly from zero, the two independent modes with angular 
momentum L = I. Computations are much easier for 
cusp or multicusp kernels than for the phenomenologi- 
cal kernel used by Aminoff et al. [ll|, because cusp ker- 
nels are easily inverted to give Green's functions. Even 
with no adjustment of the weights and sharpnesses of the 
three-cusp kernel, the modeled signal is very close to the 
observed signal shown in Fig. 4 of Aminoff et al. [3] for 
57 niTorr of Ne buffer gas. 

Fig. [5] shows the modeled prediction of what would be 
observed if the longitudinal field of Aminoff et al. [l^ 
were replaced by a transverse field of 0.1 G and with 
other experimental conditions the same as those of Fig. \5\ 
Quite different signals should be observed for magnetic 
fields large enough to cause substantial spin rotation be- 
fore the optically pumped atoms can reach the cell walls. 
This is because a spin-polarized atom generated in one 
velocity group will rotate around the transverse magnetic 
field in the time needed to diffuse (in velocity space) to 
the velocity group detected by the circularly polarized 
probe beam. It is entirely possible for the spin polariza- 
tion to rotate more than 90°, reversing the sign of the 
probe signal in the process. Magnetic fields with trans- 
verse components excite the 12 coherence modes with 
Am = ±1, and if there are also longitudinal components 
of the magnetic fields the three population modes men- 
tioned above are excited as well. So for non-longitudinal 
fields at least 15 modes \sj) and corresponding mode am- 
plitudes \xj) need to included in the sums of (|147p . and 
numerical calculations require more time. The nature 
of the "magnetic depolarization" signals generated by 
transverse magnetic fields is very sensitive to the param- 
eters (sharpness, weight) used to construct the multicusp 
kernel. Systematic experiments with magnetic depolar- 
ization of the MCD signals at various buffer gas pressures 
and various magnetic fields, analyzed with the efficient 
modeling methods we have outline in this paper, would 
provide a very good way to determine the optimum pa- 
rameters for cusp kernels. 



B. Velocity-selective optical pumping in potassium 

vapor 

In our own laboratory, we have completed preliminary 
experiments on velocity-selective optical pumping with K 
vapor. The basic experimental arrangement, depicted in 
Fig. [71 is almost identical to that of earlier experiments 
by Bloch et al. [19| which were done with no buffer gas. 
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FIG. 7. (Color online) Schematic of a saturated-absorption 
experiment, described in more detail in the text. For the 
somewhat similar experiment of Aminoff et al. , modeled 
in Fig. O and for the proposed experiment of Fig. [S] the in- 
tensity chopper is replaced by a device that modulates the 
circular polarization of the pumping light, and a weak probe 
beam has fixed circular polarization. Alternatively, one could 
modulate the polarization of the probe beam and use fixed 
circular polarization for the pump beam, which would elimi- 
nate any concerns about transient buildup and decay of the 
optically pumped spin polarization. 



The pump beam was produced by a diode laser and was 
frequency-scanned across the 770 nm Dl absorption line 
of K atoms. The probe beam was a retro-reflected frac- 
tion of the pump beam. The pump-beam intensity was 
reduced with neutral density filters to keep it well within 
the first-order regime where (|148|) describes the attenua- 
tion of the probe light. To ensure that optical pumping 
by the the probe beam was negligible, the probe inten- 
sity was attenuated by about about a factor of ten with 
respect to the pump beam. Both the pump and probe 
beams were linearly polarized along a 1 G magnetic field. 
The approximate radius of the beams was a — 1 mm. 
The pump beam intensity was modulated on and off at 
80 Hz with a chopper wheel. The pump beam produced 
a modulated spin polarization of the atoms which mod- 
ulated the transmitted intensity of the probe beam. The 
intensity modulation of the transmitted probe beam was 
detected with a lock-in amplifier, referenced to the chop- 
per wheel. 

Our glass cells included both spheres and cylinders, 
but for modeling we took a representative cylinder with 
a radius 6 = 5 mm. Before filling, a few drops of K metal 
were distilled under vacuum into the cell. The metal was 
of natural isotopic abundance, 93.26% '^^K and 6.73% 
""^K, and the relatively small signals from ''^K can be 
seen in some experimental data. Most cells were filled 
with low pressures of Kr or N2 gas but some "vacuum" 
cells had no intentional gas. Some of these vacuum cells 
were found to have small amounts of gas in them, as 
was shown clearly by collisional pedestals in the VSOP 
signals. During the experiments, the K vapor pressure 
was controlled by keeping the cell at a temperature near 
50 °C in a temperature-stabilized oven with glass win- 
dows. The optical depth of the vapor was intentionally 
kept low (about 0.1 at the peak of the Doppler profile) in 
order to keep the laser intensity relatively uniform down 
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FIG. 8. (Color online) The left panels are saturated absorp- 
tion signals from potassium vapor, measured with the appa- 
ratus of Fig. [T] The right panels are signals modeled with 
(|148|l for monoisotopic ^^K vapor. The buffer-gas pressures 
were: (a) No intentionally added buffer gas, (d) 1 mTorr of 
N2; (b) or (e) 100 mTorr of Kr, (c) or (f) 1500 mTorr of N2; 
The modeling parameters were: diffusion coefficient at 1 atm 
Do = 0.10 cm^ s~^ for both Kr and N2, laser beam diam- 
eter a — 1 mm, cell diameter 6 = 5 mm. The cusp-kernel 
parameters of Fig. [2] were used with no adjustment. 



the entire length of the cell. As shown in Fig. [SJ the dom- 
inant '^^K signals observed with the apparatus of Fig. [7] 
could be simulated very well with the first-order absorp- 
tion rate of (|148p . with the parameters mentioned in the 
figure caption. Because of the experimental arrangement, 
which cannot excite coherences, only the 8 population 
modes out of the 64 possible spin modes \sj) have ampli- 
tudes Ix;^;^) that differ appreciably from zero. 



Appendix: Spatial diffusion and velocity damping 

Here we show that the velocity damping rate 7vd can 
be inferred from the measured spatial diffusion rate D 
and the smallest, non-zero eigenvalue, ai of the collision 
operator with the formula 



7vd 



2aiD 



(A.l) 



This is consistent with Eq. (2.28) of Herman, Haverkort 
and Woerdman [20] . 

Let u be the spatial position in units of vo/jvd, 
the characteristic distance an atom can travel in one 
velocity-damping time, l/7vd- The dimensionless time 
is r = t7vdi where l/7vd is the characteristic time be- 
tween velocity-changing collisions. We ignore the atomic 
spin and consider atoms with a number density TV = 
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N{u,x,t) = ^'^^QNn{u,T){x\vn), such that at the time 
r the probabihty of finding a particle with position be- 
tween u and u + du and velocity between x and x + dx 
is N{u,x,T)dudx. If we quantize the velocity x, we can 
write the density N as an abstract column vector 



(A.2) 



with velocity- mode amplitudes — {vn\N). We as- 
sume that |A^) evolves according to the Boltzmann equa- 
tion in one dimension, 



A 



vd 



(A.3) 



The collision operator A^d was given by (j56p . We may use 
the identity xHn — nHn-i + Hn+i/2 and the definition 
(|Bn)) of the velocity modes to write the velocity operator 
X in (|X3l) as 



X 



/ n + 1 



n=0 



\Vn+l)iVn\ + \Vn)iVn+l\}. (A.4) 



Substituting (|A.4p and (|A.2p into (lA.Sp . multiplying the 
resulting equation on the left by ((wo| and ((wi|, and re- 
membering that Qfo = wc find 



dr 
dr 



1 dNi 
2~du 

Id No dN2 

2 du du 



aiNi. 



(A.5) 
(A.6) 



We are interested in "late-time" distributions N that 
have Maxwellian, or very nearly Maxwellian velocity dis- 



tributions, where the mode amplitudes iV„ = ((z;„|iV) de- 
crease rapidly with n so that TVg ^ TVi ^ iV2 ^ • • • . 
Retaining only the first two amplitudes, iVo and A^i, mul- 
tiplying (jA.sp by d/dr, multiplying (IA.6I) by d/du, and 
combining the results to eliminate iVi we find 



dr^ 



2 du^ 



ai- 



dNo 
dr 



(A.7) 



The distribution Nq will evolve more and more slowly 
with increasing time, so for sufficiently late time we can 
ignore d'^No/dr^ compared to dNo/dr in (IA.7I) and find 
the diffusion equation 



dNo 



1 d'^No 



(A., 



dr 2ai du^ ' 
with the diffusion coefficient (in dimensionless units) 



D 



1 

2q7 



(A.9) 



For a cusp kernel with sharpness s we see from (j64p that 
ai = l/{s + 1), so (|A.9P gives D ^ {s + l)/2. For large 
sharpnesses, the diffusion coefficient is very nearly half 
the value of the sharpness, which is intuitively reasonable 
since sharp kernels represent velocity damping dominated 
by grazing-incidence collisions, which do little to hinder 
spatial diffusion. Using the unit of length voHvd and 
unit of time 1 /7vd to convert the dimensionless diffusion 
coefficient (IA.9I) to dimensional units (cm^s~^ for cgs 
units) we find (jA.ll) . 
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